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Analysis of the nonlinear Schrodinger vortex reconnection is given in terms 
of coordinate-time power series. The lowest order terms in these series cor- 
respond to a solution of the linear Schrodinger equation and provide sev- 
eral interesting properties of the reconnection process, in particular the non- 
singular character of reconnections, the anti-parallel configuration of vortex 
filaments and a square-root law of approach just before/after reconnections. 
The complete infinite power series represents a fully nonlinear analytic so- 
lution in a finite volume which includes the reconnection point, and is valid 
for finite time provided the initial condition is an analytic function. These 
series solutions are free from the periodicity artifacts and discretization error 
of the direct computational approaches and they are easy to analyze using a 
computer algebra program. 
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1. INTRODUCTION 

Vortex solutions of the nonlinear Schrodinger (NLS) equation are of 
interest in nonlinear optics EEE1 anc j j n the theory of Bose-Einstein conden- 
sates^ (BEC). The NLS equation is also often used to describe turbulence 
in superfluid helium. ^ NLS is a nice model in this case because the vortex 
quantization appears naturally in this model and because its large-scale limit 
is the compressible Euler equation describing classical inviscid fluids. At 
short scales, the NLS equation allows for a "quantum uncertainty principle" 
which allows vortex reconnections without the need for a finite viscosity 
or other dissipation. Numerically, NLS vortex reconnection was studied by 
Koplik and Levine 1 ^ and, more recently, by Leadbeater et a/P and, for a 
non-local version of NLS equation, by Berloff et alW In applications to su- 
perfluid turbulence, the NLS equation was directly computed by Nore et 



S.V. Nazarenko and R.J. West 



aiP Such cryogenic turbulence consists of repeatedly reconnecting vortex 
tangles, with each reconnection event resulting in the generation of Kelvin 
waves on the vortex cores ^ and a sound emission. ™ These two small-scale 
processes are very hard to correctly compute in direct simulations of 3D NLS 
turbulence due to numerical resolution problems. A popular way to avoid 
this problem is to compute vortex tangles by a Biot-Savart method (derived 
from the Euler equation) and use a simple rule to reconnect vortex filaments 
that are closer than some critical ( "quantum" ) distance. This approach was 
pioneered by Schwarz ^ and it has been further developed by Samuels et 
o/P In this case, it is important to prescribe realistic vortex reconnection 
rules. Therefore, elementary vortex reconnection events have to be carefully 
studied and parameterized. Numerically, such a study was performed by 
Leadbeater et a/.,® the present paper is devoted to the analytical study of 
these NLS vortex reconnection events. 

The analytical approach of this paper is based on expanding a solution 
in powers of small distance from the reconnection point, and small time mea- 
sured from the reconnection moment. The idea is to exploit the fact that 
when vortex filaments are near reconnection, the nonlinearity in the NLS 
equation is small. This smallness of the nonlinearity just stems from the 
definition of vortices in NLS (curves where \P = 0) and the continuity of \E'. 
Their core size is of the order of the distance over which ^ — > 1 (where Vl/ = 1 
represents the background condensate). Therefore, for vortices near recon- 
nection, separated by a distance much smaller than their core size, VP is small 
provided it is continuous. Thus, to the first approximation the solution near 
the reconnection point can be described by a linear solution which, already 
at this level, contains some very important information about the recon- 
nection process: (1) that the reconnection proceeds smoothly without any 
singularity formation, (2) that in the immediate vicinity of the reconnection 
the vortices are strictly anti-parallel and (3) just before the reconnection 
event the distance between the vortices decreases as It] 1 / 2 , where t is the 
time measured from the reconnection moment. Note that result (1) could 
surprise those who draw their intuition from vortex collapsing events in the 
Euler equation (which are believed to be singular). On the other hand, re- 
sults (2) and (3) are remarkably similar to the numerical and theoretical 
results found for the Euler equation. I 1 ? ! 14 ! 1 ? ! 

In section II of this paper we examine the local analysis of the recon- 
nection process by deriving a linear solution and in section III consider its 
properties. The linear solution describes many, but not all the important 
properties of vortex reconnection. In particular, it cannot describe solu- 
tions outside the vortex cores and, therefore, it cannot describe the far-field 
sound radiation produced by the reconnection. On the other hand, one can 
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substitute the linear solution back into the NLS equation and find the first 
nonlinear correction to this solution. Recursively repeating this procedure, 
one can recover the fully nonlinear solution in terms of infinite coordinate 
and time series. This derivation is discussed in detail in section IV. The se- 
ries produced are a general solution to a Cauchy initial value problem. Thus, 
by Cauchy-Kowalevski theorem, ^ these series define an analytic function 
(with a finite convergence radius) provided the initial conditions are analytic. 
The generation of such a suitable initial condition is addressed in section V. 
Our series representation of the solution to the NLS equation is exact, and 
therefore will include such properties as sound emission. However, due to 
the finite radius of convergence of the analytic solution, one is unable to ob- 
serve a far-field sound emission directly. In this paper, we use Mathematica 
to compute some examples of the fully nonlinear solutions for the vortex 
reconnection. The results of which are presented in section VI. 

Let us summarise the advantages and disadvantages that our analytical 
solution has with respect to those being computed via direct numerical sim- 
ulations (DNS). Firstly, our analytical solutions are obtained as a general 
formula, simultaneously applicable for a broad class of initial vortex posi- 
tions and orientations. Secondly, our analytical solutions are not affected by 
any periodicity artifacts (which are typical in DNS using spectral methods) 
or by discretization errors. On the other hand, our analytical solutions are 
only available for a finite distance from the vortex lines (of the order of the 
vortex core size) because their defining power series have a finite radius of 
convergence. 



2. LOCAL ANALYSIS OF THE RECONNECTION 

Let us start with the defocusing NLS equation written in the non- 
dimensional form, 

i* t + A* + (l- |^| 2 )^ = 0. (1) 

Suppose that in vicinity of the point r = (x, y, z) = (0, 0, 0) at t = to we have 
Vl/ = \&o such that Re^o = z, and Im^o = az + bx 2 — cy 2 , where a, b and 
c are some positive constants. For such initial conditions the geometrical 
location of the vortex filaments, = 0, is given by two intersecting straight 
lines, z = and y = ±^b/cx. 

In the small vicinity of the point r = 0, deep inside the vortex core 
(where ~ 0), we can ignore the nonlinear term found in equation (^Q). 
Further, by a simple transformation ^ = <J>e** we can eliminate the third 
term \I/ and obtain z<J?t + A<I> = 0. (This just corresponds to multiplying 
our solution by a phase, it does not alter its properties, but does make the 
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following analysis simpler). It is easy to see that the initial condition has 
not changed under this transformation, = ^o- Advancing our system a 
small distance in time t — to, we find Re = Re \&o — (t — to) Aim and 
Im $ = Im + (t - t ) ARe <& , or 



For both t — to < and t — to > the set of vortex lines, $ = 0, is given 
by two hyperbolas. A bifurcation happens at t = to where these hyperbolas 
degenerate into the two intersecting lines (see Fig. ^). This bifurcation 
corresponds to the reconnection of the vortex filaments. Thus, we have 
constructed a local (in space and time) NLS solution corresponding to vortex 
reconnection. Obviously, this solution corresponds to a smooth function $ at 
the reconnection point. It should be stressed that this is not an assumption, 
but just the way in which we have chosen to construct our solution. However, 
we do believe that this observed smoothness is a common feature of NLS 
vortex reconnection events. If this is true then all such reconnecting vortices 
could locally be described by the presented solution as the intersection of a 
hyperbola with a moving plane provides a generic local bifurcation describing 
a reconnection in the case of smooth fields. 



3. PROPERTIES OF THE VORTEX RECONNECTION 

The local linear solution we have constructed © reveals several impor- 
tant properties of the reconnection of NLS vortices. 

1. Whatever the initial orientation of the vortex filaments, the reconnect- 
ing parts of these filaments align so that they approach each other in an 
anti-parallel configuration. Indeed, according to (j2J), the fluid velocity field 
v = Varctan([/m$]/[i?e$]) = V &rcta.n([az+bx 2 -cy 2 ]/[z-2(b-c)(t-to)]). 
At the mid-point between the two vortices one finds a velocity field consis- 
tent with an anti-parallel pair, v = l/[2a 2 (c-b)(t-t )]e z / 0. (For a parallel 
configuration one would find v = 0). Amazingly similar anti-parallel con- 
figurations have been observed in the numerical Biot-Savart simulations of 
thin vortex filaments in inviscid incompressible fluids J12HEI 

2. The reconnecting parts of the vortex filaments approach each other as 
\/i — to. Indeed, setting Re <3? = Im = and y = in (j2j) one obtains 
x = ±y/2a[(c/b) — l](t — to). Exactly the same scaling behaviour, for ap- 
proaching thin filaments in incompressible fluids, has been given by the 
theory of Siggia and Pumir l 13 l 14 l 15 l an d has been observed numerically in 
Biot-Savart computations. ECU 




(2) 
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(d) t = -0.1 (e) t = 0.0 (f) t = 0.1 



Fig. 1. Linear solution Eq. © of the nonlinear Schrodinger equation for 
a = 1, 6 = 3, c = 2 and to = 0. Sub-figures (a), (b) and (c) show the 
intersection of the real (plane) and imaginary (hyperbolic paraboloid) parts 
of Eq. © at successive times t = —0.1, t = 0.0 and t = 0.1 respectively. Sub- 
figures (d), (e) and (f) show the corresponding lines of intersection where 
\P = 0; reconnection occurs at t = to = 0. 



3. The nonlinearity plays a minor role in the late stages of vortex recon- 
nection in NLS. This is a simple manifestation of the fact that in the close 
spatio-temporal vicinity of the reconnection point ~ 0, so that the dy- 
namics are almost linear. 

This last property can be also reformulated as follows. No singularity is ob- 
served in the process of reconnection according to the solution both the 
real and imaginary parts of behave continuously in space and time. This 
property is in drastic contrast to the singularity formation found in vortex 
collapsing events described by the Euler equation. Indeed, distinct from in- 
compressible fluids, no viscous dissipation is needed for the NLS vortices to 
reconnect. Here, dispersion does the same job of breaking the topological 
constraints (related to Kelvin's circulation theorem) as viscosity does in a 
normal fluid. 
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4. NONLINEAR SOLUTION 

We will now move on to consider the full NLS equation. We will use a 
recursion relationship to compute the solution assuming that x,y,z ~ e and 
t ~ e 3 (for simplicity, we take to = 0)- The solution we obtain will therefore 

be of the form * = *(°) + *W + H , where #( n ) ~ e n . The above e 

scaling of x, y, z and t has been chosen to generate a recursion relationship 
when substituted in the NLS equation (0). Of course we could have chosen 
a different e dependence, however, as the final series representation of our 
solution contains an infinite number of terms, this would just correspond to 
the same solution but with a suitable re-ordering. 

Consider the NLS equation (JTJ. Firstly, we note that dt ~ e~ 3 and 
A ~ e" 2 and therefore i*j m) ~ e m " 3 , Af( n ) ~ e n " 2 , [|*| 2 *] (p) ~ & and 
$r(<?) ~ e q , where | a ¥]W = Eij=i Matching the terms, 

by setting m = n + 1 and p = q = n — 2, and integrating we find 

^(n+l) = ^(n+1) + ■ /" U^(n) + ^(n-2) _ [|^|2^](n-2)l (3) 
JO L J 

where ^q 1 ^ are arbitrary n th order functions of coordinate which appear as 
constants of integration with respect to time. The full nonlinear solution 

(n) 

of the Cauchy initial value problem can now be obtained by matching ^ 
to the n th order components of the initial condition at t = obtained via a 
Taylor expansion in coordinate. Let us assume that the initial condition is an 
analytic function so that it can be represented by power series in coordinates 
with a non-zero volume of convergence. Then, by the Cauchy-Kowalevski 
theorem, the function \P will remain analytic for non-zero time. In other 
words, the solution can also be represented as a power series with a non-zero 
domain of convergence in space and time. Remarkably, the recursion relation 
Eq. (j3J) is precisely the means by which one can write down the terms of 
the power-series representation of the fully nonlinear solution to the NLS 
equation, with an arbitrary analytical initial condition ^o- 



5. INITIAL CONDITION 

Our next step is to construct a suitable initial condition for our study 
of reconnecting vortices. This initial condition will have to be formulated 
in terms of a power series. We start by formulating the famous line vortex 
solution to the steady state NLS equation ^ in terms of a power series. 
Substituting * = Ae ie into Eq. (0), we find AA - A\V9\ 2 + A - A 2 = 0, 
where we have used the fact that V# • VA = and A9 = 0. We can simplify 
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this equation, since A = A(r) and therefore, AA = -d r (rd r A). However, we 
also note that |V 2 #| = -7 since d x 9 = —y/r 2 and d y 9 = —x/r 2 . Therefore, 
we have ^d r (rd r A) — A — A 3 + A = 0. We will solve this equation using 
another recursive method. We would like to get a solution of the form 
A = ao + a±r + a^r 2 + a^r 3 + • • • = J2 n A^ n \ (However, we can set ao to 
zero on physical grounds, since we require $ = at r = 0). As before 
ld r (rAi m) ) = m 2 a m r m - 2 ~ e m ~ 2 , ^ = a n r n ~ 2 ~ e™~ 2 , [A 3 ]W ~ 
and A® = a q ri ~ e 9 , where [A 3 ] GO = £? i=1 AW^'U^-^'). Again, by 
matching powers of r we can derive a recursion relationship for a n . Setting 
m = n and p = q = n — 2 we obtain 

a n = (/n-2 - a n _ 2 )/ (n 2 - 1), (4) 

where f p = [A 3 ]^/ r P. 

We should note that am = for all n. Therefore, taking a power of r 
out of our expansion for A(r) we find, 

= A(r)r w = rg(r)e w , (5) 

where g(r) = g(r 2 ) = a\ + a^r 2 + a^r 4 + • • • = Y^=i o-2n-\f 2n ~ 2 ■ Further, re 
is complex so we can write re tS = x + iy and hence our prototype solution, 
for a vortex pointing along the z-axis, is \l/ = (x + iy)g{x 2 + y 2 ). 

We can manipulate this prototype solution to get an initial condition 
for our vortex reconnection problem. Our initial condition ^0 will be made 
up of two vortices, and ^2) a distance 2d and angle 2a apart. Following 
the example of others, [Koplik et al, Ref. [8^ and [Leadbeater et al., Ref. 
we take the initial condition to be the product of ^1 and Vt^j that is ^0 = 
^1^2- One could argue that such an initial condition is rather special, as two 
vortices found in close proximity would typically have already distorted one 
another in their initial approach. Nevertheless, such a configuration provides 
us with a valuable insight into the dynamics of NLS vortex reconnections. 

Firstly, we would like the vortices in the (x, y) plane. We can do this 
by transforming our coordinates x — > y, y —* z and z — > x. This will 
give us a vortex pointing along the x-axis ^ = (y + iz)g(y 2 + z 2 ). The 
vortex can now be rotated by angle a to the x-axis in the (x, y) plane via 
x — > x cos a — y sin a and y — » y cos a + x sin a. Finally, we shift the whole 
vortex in the z direction by a distance d using z — > z — d we finally obtain 
\?i = [y cos a + x sin a + i(z — d)]g((y cos a + x sin a) 2 + (z — d) 2 ). In a similar 
manner, ^ is a vortex at angle —a and shifted by —d in the z direction, 
^2 = [ycosa — xsina + %{z + d)]g((y cos a — xsina) 2 + (z + d) 2 ). 
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Fig. 2. The initial condition for the prototype vortex solution is constructed 
via an appropriate expansion for A = A(r), Eq. Q). Here we can see the 
expansion for A = A(r) truncated at three different orders of n; n = 5 (dash- 
dot line), n = 15 (dashed line) and n = 21 (solid line) with a\ = 0.6. At 
higher orders one would see the existence of a finite radius of convergence at 
r ps 2.5. 

6. RESULTS 

It would time consuming to expand the analytical solution, derived in 
the previous section, by hand. Thankfully, we can use a computer to perform 
the necessary algebra, and to derive the hugh number of terms the recursive 
formulae generate. What follows is an example solution of the reconnection 
of two initially separated vortices. 

Firstly, we need to consider the validity and accuracy of our initial con- 
dition. Fig. [21 shows the prototype solution A = A(r) for a single vortex, 
at various different orders. Increasing the order will obviously improve ac- 
curacy. However, one should note that at higher order there is evidence 
of a finite radius of convergence r c . This will restrict the spatial region of 
validity for our full t-dependent solution. Our prototype solution also has 
a dependence on a±. In the following simulation we have chosen a\ = 0.6 
numerically so that the properties of A(r) match that of a NLS vortex. It is 
evident that we cannot satisfy these properties completely (namely — * 1 
as r — > oo) as our power series diverges near r c . Nevertheless, this does not 
present us with a problem if we restrict ourselves to considering the evolution 
of contours of such as |^| < 1, where A(r) is realistically represented. 
Further, it should be noted that sound radiation could in principle be visu- 
alized in our solution by drawing contours of |^| close to unity. However, to 
have an accurate representation, we would need to take a very large number 
of terms in the series expansion, therefore the study of sound in our model 




Analytical solution for nonlinear Schrodinger vortex reconnection 




(d) t = 0.3 (e) * = 0.4 (f) t = 0.5 



Fig. 3. Sub-figures (a) to (f) show the evolution of two initially separated 
vortices in time t. This realization is for \$>\ = 0.1, a\ = 0.6, d = 0.6 and 
q = 7r/4. The reconnection and separation events are clearly evident. 

is somewhat harder than the analysis of the vortices themselves. Of course 
the validity of the full t-dependent solution will be restricted, in the spatial 
sense, by the initial condition's region of convergence. The region of con- 
vergence will evolve, remaining finite during a finite interval of time (by the 
Cauchy-Kowalevski theorem), but then may shrink to zero afterwards. 

We will now discuss an example solution. As we only wish to demon- 
strate this method, we will not consider a high order solution in this paper. 
In our example simulation below, we used Mathematica to perform the nec- 
essary algebra in generating a nonlinear solution up to 0(e 6 )- (One should 
note that although the prototype solution © for a single vortex has a,2 n = 
for all n, our initial condition is made up of two vortices, i.e. two series 
multiplied together. Therefore, there will be cross terms of order 0(e 6 ) in 
our initial condition). 

Our choice of parameters will bed = 0.6 and a = n/A. This corresponds 
to two vortices, initially separated by a distance 1.2, at right angles to each 
other. Fig. 01 shows the evolution of the iso-surface \^>\ = 0.1 in time, 
demonstrating reconnection and then separation. Examining this solution 
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in detail we can clearly see evidence of some of the properties mentioned 
earlier - that of a smooth reconnection (the absence of singularity) and the 
anti-parallel alignment of vortices prior to reconnection. 



7. CONCLUSION 

In this paper we presented a local analysis of the NLS reconnection 
processes. We showed that many interesting properties of the reconnection 
can already be seen at the linear level of the solution. We derived a recursion 
formula Eq. (J3J) that gives the fully nonlinear solution of the initial value 
problem in a finite volume around the reconnection point for a finite period 
of time. In fact, formula (j3J) can describe a much wider class of problems. 
Of interest, for example, are solutions describing the creation or annihilation 
of NLS vortex rings. This process is easily described by considering vortex 
rings, at there creation/annihilation moment, as the the intersection of a 
plane with the minimum of a paraboloid. Further, this method of expansion 
around a reconnection point can be used for other evolution equations, e.g. 
the Ginzburg-Landau equation. These applications will be considered in 
future. We wish to thank Robert Indik, Nicholas Ercolani and Yuri Lvov for 
their many fruitful discussions. 
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